Quantum Timing and Synchronization Problems 



Diego dc Falco 
Dipartimcnto di Scicnzc dcll'Informazione 
Universita degli Studi di Milano 
via Comelico 39, 20135 Milano. Italy 
e-mail: defalco@dsi.unimi.it 



Dario Tamascclli 
Dipartimcnto di Matcmatica 
Universita degli Studi di Milano 
via Saldini 50, 20133 Milano, Italy 



e-mail: tama@mat.unimi.it 



Abstract 



Feynman's model of a quantum computer provides an example of a continuous-time 
quantum walk. Its clocking mechanism is an excitation of a basically linear chain of spins 
with occasional controlled jumps which allow for motion on a planar graph. The spreading 
of the wave packet poses limitations on the probability of ever completing the s elementary 
steps of a computation: an additional amount of storage space 5 is needed in order to 
achieve an assigned completion probability. In this note we study the END instruction, 
viewed as a measurement of the position of the clocking excitation: a 7r-pulse indefinitely 
freezes the contents of the input/output register, with a probability depending only on 
the ratio S/s. 

Keywords: continuous-time quantum walk; quantum END problem; telomeric chain; 
7r-pulse trap; quantum subroutines; Grover's algorithm. 

1 Introduction 

It has been shown by Feynmanfl] that it is possible to implement the sequential application, in 

the desired order, of a collection A = A s _i A2 ■ A\ of unitary operators to an input/output 

register by using s additional program counter sites. 

For the sake of definiteness, we will think of each program counter site j = 1,2,... ,s as 
occupied by a spin-1/2 system r(j) = (ri(j), T2(j), T$(j)). We will refer to the collection of 
such spins, which act in effect as a quantum clocking mechanism, as to a "program line" . 
The input/output register will be, similarly, implemented by a collection of a certain number 
/i of spin-1/2 systems a(i) = (ci(i), <72(i), 03(2)), i = 1,2,...,//. 

The evolution of the complete system, register + program line, will be given by the Schrodinger 
equation: 



where the Hamiltonian H is supposed to be time independent and involving at most 
3-body interactions: 
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We have indicated by t + (j) and T—(J) the raising and lowering operators for the third 
component of the spin occupying the j-th program counter site. 

Since it is requested that each term of the sum is (at most) a 3-body interaction, each Aa 
either acts on a single spin of the register or is a constant: as shown in Ref.pP, this does not 
restrict the class of functions computable by the model. 

We will restrict our considerations to initial conditions of the form: 

\ip(0)) = {register^) <g> \program lineo), (3) 

where {program lineo) — | T3 ( 1 ) = +l,T%{j) = —1 for 1 < j < s) describes the state in 
which only the spin located at program site number 1 is "up" . 

It helps the intuition to think of the initial state \register$) of the register as a simultaneous 
eigenstate of the components of the a spins in selected directions, encoding the initial word 
(or superposition of words) on which the machine is required to act. 
Because of the conservation law [H, N3] = 0, where 

N .^±i±^il, ( 4) 

the above initial condition has the important consequence that the evolution \ip(t)) takes 
place in the 2^s dimensional eigenspace of N3 belonging to the eigenvalue +1. 
The intuition of "a single clocking excitation travelling along the program line" emerging 
from the above conservation law is made precise by introducing the observable position of the 
excitation, or position of the cursor: 

Q't^-^f 1 - (5) 

5=1 

It is then easy to convince oneself that the evolution of the overall system is of the form 

s 

\tp{t)) = c(t, k; s) \register k _i) ® \Q = k) (6) 

k=l 

where, for 1 < h < s — 1, \registerh) = Ah\registerh~\) ■ 
In Feynman's words (adapted to our notations), (jEJ) says that, starting from the initial con- 
dition "If at some later time the final site s is found to be in the \ts(s) = +1) state 
(and therefore all the others in \T^(j) = — 1) ), then the register state has been multiplied by 

A s -i A2 ■ A\ as desired". It as been shown in Ref.j2] that this is a somewhat big "If", 

under two respects: 

i. at no instant of time the probability |c(t, s; s)\ 2 is larger than const ■ s~3 ; 

ii. the cursor keeps bouncing back and forth between positions 1 and s, thus in effect 
making the above upper bound attainable only at selected instants of time. 

The above two statements are reviewed and made quantitative in Section |2j 

Section |31 is devoted to the "quantum END problem" : we remove the cursor in order to 

prevent it from returning down the program line and "undoing" the computation. Removing 
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the cursor and storing the result of the computation in the contents of the register is in effect 
a measurement procedure, that in Sec. |3] will be modelled by a suitable time dependent 
perturbation (a 7r-pulse) applied to a variant of the Hamiltonian 
Section |I] is devoted to conclusions and outlook. 

2 The Motion of the Cursor 

We recall, first of all, that the motion of the cursor does not depend on the operators acting 
on the register, pp For the particular case of a sequential program line as the one described 
by the Hamilonian ©, this is made evident by the explicit expression [3] of the amplitudes 
c(t,k;s) in Q: they are given, independently of the operators Aj, by: 

2 s 

c(t,k;s) = } exp (it cos (0(n; s))) sin (6>(n; s)) sin (k 9(n; s)), (7) 

s + 1 ^— ' 

n=l 

where: 

s nit , . 

0(n;s) = — 8 
s + 1 

Similar results hold in the case, studied in Ref.jl], in which, because of conditional jumps 
in the program line (such as the ones needed in the iteration of quantum subroutines), the 
cursor performs, in effect, a continuous-time quantum walkjH] on a planar graph. In this note 
we restrict ourselves to the sequential case. 

The main purpose of this Section is to give examples of the behaviour recalled in the obser- 
vations (0) and (juj) of SecO 

This we do with the help of the following Hamiltonian: 

s-l 

H = --(XV+(j + l)A 7 -r_(i) + 
i=i 

+r+(s+ l)/9_T_(s) + 
s+5-l 

+ E 7-+(j + l)r_(j) + ^c). (9) 
i=s+i 

With respect to the Hamiltonian (J2J), we have introduced an additional control q-bit p = 
(pi, P2, Ps) in the term r+(s + l)p_r_(s); this is an example of a conditional jump in the 
quantum walk performed by the cursor: it acts non trivially only in the eigenspace belonging 
to the eigenvalue +1 of ps, enabling the transition \Q = s) — > \Q = s + 1). If this transition 
is enabled, then the cursor can visit the additional telomeric sites s + 1, . . . , s + 5, else it gets 
reflected back. 

Figures n and |21 give examples of the behaviour of the probability 

s+S 

P (s<Q) (*) = P(,<Q<s+S)(t) = E W*>-7')! 2 ( 10 ) 

3=8 

of finding the register in the state A = A s _\ A<i ■ A\\register^) , under two different 

initial conditions, which determine two different forms of the amplitudes 7. 
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Figure ^ corresponds to the initial condition \program line®) = \p% = — 1)<8>|Q = 1): the mo- 
tion of the cursor remains confined to the sites 1, . . . , s, as it is j(t, k) = c(t, k; s) if 1 < k < s, 
otherwise. The probability P( s <Q)(t) of finding the computation completed satisfies in this 
case the inequality: 

P(.<Q)(t) < 4 (11) 



S 3 



Figure [2 corresponds to the initial condition: \program lineo) = \ps = +1)<8> \Q = 1), leading 
to j(t, k) = c(t, k; s + 5) for 1 < k < s + <5. For t just below s + 25 the probability P( s< Q\(t) 
of finding the computation completed is close to the much less severe upper bound: j^j 



P (s<Q<s+8)(t) < 1 



2 I ■ ( 1 
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Figure 1: \program lineo) = \P3 = — 1) <8> |Q = 1); s = 20. 
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Figure 2: {program lineo) = \P3 = +1} Cg) |Q = 1); s = 20; 5 = 10. 



3 The Quantum END Instruction 

The abrupt collapse of P( S <Q) (*) a ^ time t « s + 25, evident from Fig. |21 corresponds to the 
following fact: travelling with average speed close to 1, at time t ~ s + 25 the cursor "returns 
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down the active part of the program line", thus, in effect, undoing the calculation. 
Bringing the computation to an END, and storing the result is not completely trivial in 
the case examined here of a reversible quantum clocking mechanism: "Surely a computer 
has eventually to be in interaction with the external world , both for putting data in and for 
taking it outQ]". 

A simple model of such interaction is suggested by inspection of the Hamiltonian @ : starting 
from the initial condition \p$ = +1), the transition \Q = s) — > \Q = s + 1) is enabled by the 
control term r+(s + l)p_r_(s) which, simultaneously, determines the transition \p% = +1) — > 
\P3 = "I)- 

The transition \Q = s+1) — > \Q = s), enabled by the hermitian conjugate term T + (s)p + T-(s+ 
1), will be therefore inhibited if, by external means, we enforce the transition \p$ = —1) — > 
|p3 = +1) at a time, close to to = s + 25, when most of the probability mass is in the region 
s,. . . ,s + 5. 

Figures 01 (where Figs. ^ and [2] are also reproduced for comparison purpose) presents the 
effect of the addition to © of the time dependent perturbation 

h(t)= l -B{t)- Pl (13) 

where the "magnetic field" B(t) is non vanishing only in a unit time interval around to, in 
which it takes the value n. With a probability depending only on the ratio S/s (see (|12|l) 




Figure 3: The solid line represents the probability of finding the cursor in the telomeric chain 
using a 7r-pulse applied at time to = s + 25. The dashed lines correspond to Figs^and[21 

between the lengths of the active part and the telomeric part of the program line, the 7r-pulse 
(|13|) definitively prevents the cursor from undoing the computation. 

The idea of a ir-pulse trap just presented works only if the control q-bit is initialised in the 
|p3 = +1) state. It is immediate to convince oneself that the following double trap Hamiltonian 
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does not suffer from the above limitation: 
H = h(t) + 

s-l s+5-1 s+25-1 

-2(^ T+(i + 1) ^' r - (j)+ T+(j + l)T-(j)+ T+(j + l)T-(j) + 

j=l j=s+l j=s+S+l 

+t + (s + l)p_r_(a) + r+(s + 5 + l)p+r_(s) + h.c)j . (14) 

With any initial condition for the control q-bit, under the action of the above Hamiltonian, 
the |p3 = +1) component of the state gets definitively trapped in the first telomeric region 
{s + 1, . . . , s + 6}, the |/93 = —1} component in the second one {s + S + 1, . . . , s + 25}. 

As a final remark of this section, we observe that, acting in effect as a Stern-Gerlach apparatus 
providing space separation between two different spin states, the term 

switch = (t+(s + l)p_r_(s) + t + (s + 5 + l)p+r_(s)) + h.c. (15) 

can be used also to model the preparation ( "putting the data in" ) of a register q-bit in a 
given spin state. 

4 Conclusions and Outlook 

Feynman's time honoured model of a quantum computer and its modern streamlined version, 
the continuous-time quantum walkj^j paradigm (in which the quantum motion of the cursor 
on a graph is the computation, irrespective of any action on the register), provide a fascinat- 
ing physical context in which to think of "time" under a quantum perspective and are rich 
sources of open problems. 

Inequality 1)11(1 is, for istance, a simple consequence of the spreading of the wave packet in the 
inertial motion (J7J) of the cursor on a finite lattice. Is inequality 1)11(1 strictly dependent on 
the model adopted here or is it a hint of a model independent "probability vs. computational 
complexity" uncertainty relation? 

Inequality 1)12(1 and the related discussion of Sec. 01 set bounds on the minimun amount of 
additional space needed in order to have a preassigned probability of storing the result of the 
complete calculation. Similarly, the preparation of a given input requires time and space in 
order to perform the required "Stern-Gerlach" preparation. Is it realistic to neglect this costs 
as it is done in the conventional performance analysis of performance of quantum algorithms? 

As a concluding remark, we wish to point out a case, Grover's algorithm[H], in which the 
nature, classical or quantum, of the clocking agent that successively applies the required 
primitives does make a difference. 

Grover's algorithm poses the problem of estimating the parameters a = (ai,...,a^) € 
{— lj+l}^ appearing in the transformation A = 1 — 2P & that an oracle is able to apply 
to the register. Here we have indicated by P a the projector on the state |o"3(i) = flj, i = 
1, ...,n) = |a 3 ). 

The estimation procedure of a starts from \registero) = \o~i(i) = +1, i = 1, . . . = |li) and 
proceeds by alternating the action of A with the action of B = 1 — 2|li)(li|. 
It can be shown [Jj that: 

|(a 3 |(B • A) n |li}| 2 = sin 2 ((2n + 1)0), with = arcsin^/ 2 ). (16) 
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The sharp maximum of (|16|) at n — n op ti ma i ~ 1 2^/ 2 keeps reappearing periodically if the 
computation proceeds indefinitely after n optima i steps. 

In Refs.jH] and the oscillatory nature of the overlap (|16ft between the target state |a3) and 
the current state (B • A) n |li) has been nicely explained in terms of the spectral gap of an 
analogue Hamiltonian of the form — 7^3) (as| — acting on the register viewed as an 

isolated system. 

In the context of the model considered in the previous sections, Grover's algorithm corresponds 
to the execution of a program line of the form ((2j) with Aj = A for odd j, Aj = B for even j. 
Starting form the initial condition |^(0)) = |li) (g> \Q = 1), and taking, for the sake of 
definiteness, an odd value of s, s = 2g + 1, it is easy to check that the overlap probability is, 
in this case, given by 



(^(t)|P a |V(t)) = ^sin 2 ((2n + l)^)(| C (t,2n + l; S )| 2 + |c(t,2n + 2; S )| 2 ) = 

n=0 

s 

= ^2\c{t, x; s)\ 2 sin 2 ^ x odd ), (17) 

x=l 

where x oc id is the largest odd number not larger than x. 
Direct inspection of (|17j) shows that (|16j) gives the conditional probability of finding the 
register in the target state, given that the cursor is in the state \Q = 2n + 1); in loose terms: 
the oscillatory behaviour (|16[) describes the overlap probability as a function of the "machine 
time", namely the position Q of the clocking excitation. 
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Figure 4: fi = 5; s = 2^ +1 + 1 = 65. 

The right hand side of ()17j) gives, instead, the overlap probability as a function of the 
"Galileian time" t, the independent variable in 

Figure 0] shows the behaviour of (^(^l-Pal^i)) as a function of t: the coupling of the register 
with the many clocking degrees of freedom has a damping effect on the oscillation of the state 
of the register. 

In the example just discussed we have used, in fact, a number of clocking degrees of freedom 
growing exponentially with fi. It is shown in Refs.jl] and that the above pseudo-dissipative 
effect can still be observed if, by a careful use of quantum subroutines, the number of program 
line sites is reduced to a polynomial in \x. 
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